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ABSTRACT 

We re-examine the correlation between the colors and the inclinations of the Classical Kuiper 
Belt Objects (CKBOs) with an enlarged sample of optical measurements. The correlation is strong 
(p = —0.7) and highly significant (> 8er) in the range 0° — 34°. Nonetheless, the optical colors are 
independent of inclination below as 12°, showing no evidence for a break at the reported boundary 
between the so-called dynamically "hot" and "cold" populations near w 5°. The commonly accepted 
parity between the dynamically cold CKBOs and the red CKBOs is observationally unsubstantiated, 
since the group of red CKBOs extends to higher inclinations. Our data suggest, however, the existence 
of a different color break. We find that the functional form of the color-inclination relation is most 
satisfactorily described by a non-linear and stepwise behavior with a color break at « 12°. Objects 
with inclinations ^12° show bluish colors which are either weakly correlated with inclination or are 
simply homogeneously blue, whereas objects with inclinations < 12° are homogeneously red. 
Subject headings: Kuiper Belt - methods: data analysis - solar system: general 



1. INTRODUCTION 

The Kuiper Belt is a disk of icy bodies having semi- 
major axes larger than that of Neptune. Its members are 
usually known as Kuiper Belt Objects (KBOs) or Trans- 
Neptunian Objects (TNOs). The distribution of their or- 
bits is structured leading to the identification of several 
dynamical families. Resonant KBOs are those which are 
trapped in mean-motion resonances with Neptune (those 
trapped in the 3:2 resonance are also known as Plutinos). 
Scattered KBOs, also known as Scattered Disk Objects, 
are essentially highly eccentric KBOs under strong gravi- 
tational influence of Neptune. Classical KBOs (CKBOs) 
possess relatively circular orbits that arc neither located 
in any strong mean-motion resonance with Neptune nor 
strongly subject to its gravitational influence. 

Since their discovery in 1992 more than 1200 KBOs 
have been identified. Due to their faintness only about 50 
can be spectroscopically studied with the currently avail- 
able instruments. Multicolor photometry provides, how- 
ever, a first-order approximation of their spectra, hence 
of their surface composition. Most KBOs can be studied 
photometrically and about 230 objects have at least one 
measured color. Their surface colors have shown to be 
most diverse, ranging from neutral and even slightly blue 
(relative to the Sun) to extremely red, s uggesting a large 



compositio nal diversity (see review by Doressoundiram 
eTaIp008 |. 

The origin of the color diversity remains unclear. Var- 
ious suggestions h ave been made in the context of colli 



sional resurfacing (Luu fe Jewitt||1996 
Delsanti et al.|2004lf ~ 



Gil-Hutton||2002 

Nevertheless, none of the proposec 
models has been able to consistently explain the colors 



Jewitt fe Luul|2001 



elsanti et al.|| 2004p . 



Thebault fe Doressoundiram| 2003 
Another possibility is that the ob 
served color differences reflect primordial compositional 
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variations (e.g. |Tegler et aL 2003). Such compositional 
differences would be hard to explain if KBOs formed in 
situ, since the temperature difference between 30 and 50 
AU is a very modest ~ 10 K. However, larger temper- 
ature and compositional differences might be possible if 
the KBO population, or part of it, did not form in place. 
Some dynamical mo dels, in fact, suggest outward m igra- 



tion of KBOs (e.g. Malhotra 1995 Gomes 20031. Al- 



though there are no detailed chemical studies to address 
how varied such compositions would be and how these 
would reflect in the surface colors, these dynamical mod- 
els imply a link between the current orbital inclinations 
of classical KBOs and their presumed location of origin. 

As a whole, KBOs do not show significant correlations 
between their colors and orbital parameters such as semi- 
major axis or perihelion distance. On the other hand, the 
CKBOs do show a correlation between orbital inclination 
and op tical color ( Tegler fc Romanishin|2000 Trujillo & 



Brown 



2002 1 and a correlation between perihelion and 



color ( [Tegler fc Roman ishin 2000 Peixinho et al. 2004 1 



In parallel, several works have pointed out the existence 
of two groups of CKBOs with a separation at ss 5° in in- 
clination. The two groups, usually referred to as "cold" 
(* 5°) and "hot" (i > 5°), have been identified us- 



ing n ot only orbital properties ( Brown|2001 Elliot et al. 
2005|) but also physical properties such as size and bina- 
ntylLevison fc Sternf200lf |Peixinho et al. [20041 |Gulbis 
et al.||2006[ |Moll et al.||2008p . Given its potential impor- 
tance, we re-examine the color-inclination relation using 
a new data set and appropriate statistical tests. 

2. DATA SET 



We use the system of Lykawka & Mukai ( 2007 1 to se- 
lect CKBOs for our sample. These authors classify the 
KBO families based on 4 Gyr dynamical simulations. 
The Classical KBOs have semi-major axes in the range 
37 < a < 48 AU and perihelion distances q > 37 AU, 
and must not be not located in any strong resonance 



2 
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(3:2, 5:3, 7:4, and 2:1). All orbital elements were gath- 
ered from the Minor Planet Center 1 . We have computed 
the orbital inclinations relative to the Kuiper Belt Plane 
(KBP ), hereafter denoted by ik, as defined by Elliot et al. 
( 2005 ) . We conducted our statistical tests using both the 
raw inclinations and ik, finding no significant difference 
between them. In the remainder of this work, we present 
all results in terms of ik- 

For our data set we have gathered all the B — R colors 
of CKBOs available in the literature or online. Several 
objects have colors reported in more than one work but 
the different measurements have been shown to be es- 
sentially compatible. Hence we chose to take the CKBO 
colors sequentially from the works that have been pre- 
senting lower and less dispersed error bars to those with 
(slightly) larger and more dispersed error bars. There- 
fore, firstly we gathered the CKBO colors from Tegler, 
Romanishin and Consolmagno's data sample 2 (Tegler 



et al. 



the 



2003") and references therein); secondly those from 
lISO Large P rogram on Centaurs and TNOs" ( Peix- 



inho et al.||2004[ and references therein); thirdl y those 
from the "Meudon M ulticolor Survey (2MS)" ( |Dores- 
soundiram et al.||2005[ and refe rences therein); fourthly 
those fro m | Jewitt et al.| (|2007[); an d lastly those from 
the online |Hamaut fe Delsairflfffl>002| MBOSS database 3 . 
The resulting CKBO sample used here has colors in the 
range 0.99 ^ B — R ^ 1.94. For reference, the col or of 
the Sun is (B - R) e = 0.99 ( |Hartmann et al.||1990 l. 

A histogram of the B — R error bars of the gathered 
CKBOs shows a rather continuous but very skewed dis- 
tribution, from 0.01 up to 0.21 peaking around 0.06 and 
with a mean value of 0.09. Six wayward objects pos- 
sess errors between 0.28 and 0.40, though, and we chose 
to eliminate them. The subsequent data consist of the 
B - R colors of 71 CKBOs (see Table [T}. As discussed 
in the next section, two data points appear to be out- 
liers: 2001QY 297 and 1998WV 24 . We chose to discard 
both and the final data set under analysis consists of 
the B — R colors of 69 CKBOs. The effects of keeping 
these two wayward objects in the sample are discussed 
in Section l3~6l 



3. DATA ANALYSIS 

A visual inspection of the B — R colors of CKBOs ver- 
sus ik in our data set instantly shows a trend between 
these two variables (see Fig.[T]). Using the statistical tools 
implemented in IDL we have analyzed this trend quan- 
titatively. The Spearman-rank correlation coeffi cient, p, 
for the total of N=69 CKBOs ( |Spearman||"l904| ) is: 



n - -0 7f)+ - 09 
p — \j.iu_ 007 



SL>8a 



(1) 



where SL is the significance level in standard deviations 
of a Gaussian probability distribution — error bars are 
estimated from 1000 bo otstrap extractions correc ted for 
non-Gaussian behavior ( |Efron fe Tibshirani||1993[ ). This 
is a highly significant correlation consistent with the pub- 
lished values. The square of the correlation coefficient, 
usually called the "coefficient of determination", gives 
approximately the proportion of the variation in the de- 
pendent variable that can be predicted by the changes 

1 http:/ /cfa-www. harvard.edu/iau/lists/TNOs. html 

2 http:/ /www. physics. nau.edu/~tegler/research/survey.htm 

3 http:/ /www. sc. eso.org/~ohainaut/MBOSS/ 
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Fig. 1. — B — R colors versus inclination to the Kuiper Belt 
Plane, ij.[deg], of our data set of 69 CKBOs (black dots). The 
two apparent outliers that were discarded from our analysis are 
indicated (empty circles). 



in the values of the independent variable. So, from 
p 2 = (— 0.70) 2 = 0.49 we may say that about half of the 
color variability can be accounted for by differences in 
orbital inclination. The other half is color variability un- 
accounted for by inclination differences and presumably 
related to some other undetermined variable or effect. 

Fig. [T] suggests that the color-inclination trend might 
be not linear: the colors of low inclination objects do not 
seem to correlate with inclination. When dividing the 
data set in two groups in inclination with equal number 
of objects we have a low inclination group with 34 objects 
(ik < 5°) and a high inclination group with 35 (ik ^ 5°). 
While the high inclination group still shows a strong and 

h0.05 
-0.04 



significant color-inclination correlation: p 



-0.81 



(SL > 8er), the low inclination one does not show any 
significant correlation: p = — O.IS^q'^J (SL = 0.7 a). 

In Fig. [2] we have drawn histograms of the B — R col- 
ors for the 34 objects with ik < 5°, for the 46 objects 
with ik < 12°, and for all objects. From this figure 
it seems that the color distributions for ik < 5° and 
ik < 12° are the same while only for ik ^ 12° do we 
start to see a significant number of blue objects. The 
color differences between two grou ps of objects ma y be 
analyzed using the Wilcoxon Test ( |Wilcoxon|[l945| , the 
non-parametric equivalent of the t-'l'est, also known as 
Wilcoxon Rank- Sum Test. The test ranks the full set 
of colors and assesses for incompatibility by comparing 
the ranks assigned to the members of each group. Com- 
paring the 34 objects with ik < 5° and the 12 objects 
with 5 ^ ik < 12° shows no evidence for color differences 
between the two groups (the significance level of incom- 
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Fig. 2. — Histograms of B — R colors with i k < 5° (top), with 
ik < 12° (middle), and all i k values (bottom). Only for ^ 12° 
do we see a significant number of blue objects, whereas objects 
between 5° and 12° do not appear different from those with i k < 
5°. 

patibility is 0.8 a). On the other hand, comparing the 46 
objects with if. < 12° and the 23 objects with i k ^ 12° 
shows a color incompatibility at a 6.3 a significance level. 
Next we study: (i) how the correlation coefficient varies 
with the inclusion of more highly inclined objects, (ii) 
how the mean colors vary with inclination, and (iii) which 
functional form best describes their behavior. 

3.1. Correlation as a function of inclination 

To further investigate the variations of the color- 
inclination trend we have successively computed p for 
CKBOs below a critical inclination cutoff i c k varying from 
3° in increments of 0.5° up to 20°. These two extrema 
were imposed so as not to calculate correlation values 
for very small sub-samples which were already out of the 
region of interest. Table [2] lists the results for each incli- 
nation cutoff i k , both for objects with inclinations below 
i k and those above i k — error bars and significance levels 



are also indicated. We see that p varies rather erratically 
until i c k = 12°, increases systematically with the inclusion 
of objects above that point, and reaches the 2 a (95%) 
typical minimum statistical threshold to have "reason- 
ably strong evidence" for correlation at i k = 13.5°. Such 
behavior suggests the presence of a homogenous set of 
colors below ik ~ 12° — 13.5° consistent with Fig. [2] 
While for the 46 CKBOs with i k < 12° we have no 
apparent color-inclination correlation (p = —0.15~t®'\®, 
SL = 1.0 a), for the 23 objects with i k ^ 12° a signifi- 
cant correlation is present (p = — 0.62+q n, SL = 3.2 a). 
We note that after moving the critical inclination cutoff 
by just 0.5° the 21 CKBOs with i k ^ 12.5° no longer 
show the canonical 3 a level correlation (p = —0.55^^4, 
SL = 2.6a), and for the 17 objects with i k ^ 14.5° 
the significance level drops below 2a (p = — 0.451q'2i)- 
Thus, the data provide no formally significant evidence 
for correlation among objects with i k ^ 14.5°. 

This first analysis shows that the CKBOs of smallest 
inclination are homogeneous and red, as other works have 
reported, but that homogeneity extends at least up to 
i k « 12° - 13.5°, not only up to i k w 5°. Further, the 
rapid decrease in the correlation found by removing the 
objects between 12° and 14° may suggest two separate 
groups of objects, each one having no color-inclination 
correlation whatsoever, populating two distinct parts of 
the Classical Kuipcr Belt. We will address this possibility 
further ahead. 

We are aware, though, that when considering objects 
with if. < i k and i k ^ i k separately we are also reducing 
the inclination spans under analysis, i.e., constraining 
the range of inclination values. We saw previously that 
only about half of the color variability can be explained 
by inclination differences. Consequently, the weakening 
of correlation values and their significance levels seen 
when splitting the data set in two inclination groups, 
could simply be a consequence of using an inclination 
range too narrow to detect any significant trend. To in- 
vestigate this possibility we analyze how the mean colors 
of CKBOs vary with inclination. 

3.2. Color differences as a function of inclination 

Evidently, if the sample shows a color-inclination trend 
the mean colors of CKBOs with i k < i k must be different 
from those with if. i k . That is, they must be statisti- 
cally incompatible. If the trend was approximately lin- 
ear, evidence for color incompatibility would simply vary 
smoothly with the number of objects above and below i k , 
as it also depends on that number. However, if there is a 
homogenous group of colors below some i k value then a 
maximum of color incompatibility between objects above 
and below that i k is expected to occur. 

Using the Wilcoxon Test, we successively compare the 
mean colors of CKBOs having i k < i k with those having 
ik ^ i%, varying i k from 3° to 20° in increments of 0.5°. 
Results for each i k are listed in the last column of Table 
[2] — the Wsl value is the significance level in standard 
deviations of a Gaussian probability distribution. The 
mean values are also indicated. The significance of these 
differences peaks at i k = 12.0°, with a value of 6.3 a. 
These results corroborate the existence of a homogenous 
set of colors below if. « 12°, as suggested by the analysis 
in the previous sections. 
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Fig. 3. — Histogram of the probability density distribution of best 
critical inclination cutoff if for maximum color differences (from 
Wilcoxon Tests), obtained by bootstrapping our data sample. The 
maximum color difference is found at if = 12.0° 
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3.3. Confidence intervals for critical inclination cutoff 

The finding of the critical inclination cutoff i k = 12° 
that separates the red group of CKBOs from the more 
blue ones, carried out in the previous sections, assumes 
that our data set is a representative sample of the CK- 
BOs. As with the correlation coefficients case, we may 
use bootstraps to estimate the confidence interval (error 
bar) of the best inclination cutoff i k obtained from the 
Wilcoxon Tests. We have made 1000 bootstrap extrac- 
tions from the data set and for each extraction we have 
looked for the i c k values that maximized the color dif- 
ferences between objects above and below it, as done in 
the previous section. Since in our analysis the i k is not 
continuous but discrete (with 0.5° steps) the bootstrap 
distribution is likely to be jagged. To avoid jaggedness a 
Gaussian noise with a — 0.25° was added to the inclina- 
tions of each bootstrap extraction (smooth bootstrap) . 

The probability density distribution of best critical in- 
clination cutoff ii for maximum color differences (from 
Wilcoxon Tests) is shown in Fig. [3] The i k that maxi- 
mizes color differences is well centered around 12°. Its 1 a 
confidence interval (68.3% percentile) is i k = 12. !^. 
The probability density distribution is not smoothly bell- 
shaped and two other small solution spikes arc also 



present: 5.8% probability for 



7.5 



o+O.O 
-0.5 



and 9.7% 



probability at i k — 14.5°lg'5. However, the associated 
probabilities of these spikes are low and they do not war- 
rant further attention. 

3.4. The color-inclination relation 

Having established that the CKBOs with ik < 12° con- 
stitute a group of homogeneously red objects, we next ex- 
amine the variation of B— R at larger inclinations and the 
apparent stepwise behavior at the edge of the homoge- 
neously red group. We consider three different functional 
forms for B — R color as a function of inclination: a) lin- 
ear; b) two-constant stepwise; c) constant-linear stepwise 
(see Fig. [I}. 



3.4.1. Linear fit 
Firstly, we performed a simple linear fit to the data, 



(B - R) = mi k + (B — R) 



(2) 



where m is the linear slope and (B — R) is the in- 
tercept. We have used a non-weighted Levenberg- 
Marquardt least-squares fit ( Levenberg|1944| Marquardt 
19631. We chose not to weight the data points using 
their error bars as those refer to the precision of each 
color measurement and not to the expected departure 
from the global trend. We have obtained the solution: 



with a x 2 
Fig. [4]). 



(B-R) = -0.0182 i k + 1.774 (3) 
1.208 and df = 67 degrees-of-freedom (see 



3.4.2. 



Two-constant stepwise fit 

Secondly, since our previous analysis also showed the 
possibility that CKBOs consist of two different homoge- 
nous groups of objects, neither with any color-inclination 
correlation, we fitted the data with a two-constant step- 
wise function: 



(B - R) = (B - R) ol ^ i k < i% 
(B-R) = (B- R) o2 ^i k >i% 



(4) 



i.e., a stepwise function with a constant color value below 
a given critical inclination i k , and with another constant 
value above i k . We have fitted this function to the data 
iteratively, changing i k from 3° to 20° with increments 
of 0.5°. For each iteration, (B — R) a i and (B — R) 2 are 
fitted while the critical i k is kept fixed. Table |3]shows the 
results obtained for each i k value. The best fifTaefined as 
the one which minimizes \ 2 , is obtained when i k — 13.0°: 



(B-R) = 1.701 <J= i k < 13.0° 
(B - R) = 1.316^4 > 13.0° 

with x 2 = 1-349 and df = 66 (see Fig. [Z}. 



(5) 



3.4.3. Constant-linear stepwise fit 
Thirdly, we chose to fit a constant-linear stepwise func- 



tion: 



(B-R) = (B- R) ol 

(B - R) = mi k + (B - R) o2 - 



ik < i c k 
ik > i c k 



(G) 



where we have a (B 
critical inclination i 



k- 



- R) i constant value below some 
and a linear behavior with slope 
m and intercept (B — R) 2 above i k . As for the previous 
case, we have fitted this function iteratively changing i k 
from 3° to 20° with increments of 0.5°. For each iteration 
(B — R)oi, m, and (B — R) 2 are fitted while i k is kept 
fixed. The fitting results for each i k are shown in Table 
[3| The best fit, from the minimum x 2 , is obtained when 



*fe = 12-0° 



(B-R) = 1.712 

(B - R) = -0.0159 i k 



1.703- 



i k < 12.0° 
ik > 12.0° 



(7) 
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Fig. 4. — Three functional forms studied as possible color-inclination behaviors of CKBOs. Left: Eq.pj] — the linear fit; center: Eq. [5] — 
the two-constant stepwise fit with i k = f3°; right: Eq.[7] — the constant-linear stepwise fit with i k = A \ 2 analysis shows that Eq.[7] 
is the best fit. 



with x 2 = 1-100 and df = 65 (see Fig. E}. The smallest 
X 2 value points to this solution as the Dest. Note that 
in our analysis the x 2 does not tell us the probability 
of eventually obtaining a better fit if we had another 
sample of CKBOs (with smaller error bars, for example) . 
This last solution is the best relative to the other cases, 
and validates our findings as discussed in the previous 
sections. 

3.5. Sharp boundary between groups or overlap? 

We have seen that CKBOs up to ik — 12° are homo- 
geneous in terms of their B — R color (see Fig. [21 and 
that they are redder than CKBOs with ik > 12°. If 
this is due to the existence of two independent popula- 
tions, then we might expect some mixing of the two at 
inclinations close to the boundary. Such mixing might 
also be expected from dynamical considerations as cold 
(low-ife) CKBOs may be pumped to higher inclinations 
(10° ~ 15°) due to interactions wi th resonances or even 
with a potential "planetoid" ( e.g., Kuchner et al 
Lykawka & Mukai|2007[ [20081) 



2002 



We cannot compute the x from the superposition of 
two functions for direct comparison with the fits obtained 
in the previous section. However, we have used the func- 
tional core-halo i nclination decomposition proposed by 
Elliot et al. ( 2005 1 to consider the possibility that the ob- 
served color systematics (Fig. [T|) result from the overlap 
of a red core population (low inclinations) with a blue 
halo population (high inclinations). Simulations show 
that this solution fails in the sense that it cannot repro- 
duce the color jump observed at 12°. The combination 
of the broad i nclination distribu tion for the halo objects 
(see Fig. 20 of |Elliot et al.||2005 1 and the large color dis- 
persion we observe for the bluer objects (la = 0.20; see 
also Fig. [2| results in a very smooth and broad color dis- 
tribution at all inclinations, except for ik < 3° where red 
objects are slightly more abundant. 

Interestingly, the distribution of B—R colors of CKBOs 
below ik = 12° is remarkably Gaussian, which attests to 



their color ho mogeneity. The Ko l mogorov-Smirn ov Test 
(hereafter KS; Kolmogorov 1933 Smirnov ] [l939l ) gives a 
confidence level of 99.3% (2.7 a) that the colors of ik < 
12° CKBOs are drawn from a Gaussian distribution with 
mean /i = 1.71 and standard deviation a = 0.11 (values 
calculated from the sample). This stands in contrast to 
the color distribution of objects above ik — 12° (mean 
[i = 1.34 and standard deviation a = 0.20) for which the 
KS Test gives a probability of only 22.5% of being derived 
from a Gaussian. From the KS Test the whole sample of 
B — R colors of CKBOs has only a 11.4% probability of 
being Gaussian. 

The lack of a break in the color distribution near 5° 
stands in sharp contrast to the reportedly bimodal distri 



2001| Elliot et al. 



bution of orbital inclinations (Brown 
|2005[ ). Some models have attempted to relate dynam- 
ically cold (i < 5°), red KBOs to a primordial trans- 
Neptunian disk source while dynamically hotter (i > 5°), 
blue KBOs are supposed to originate by outward scat- 
tering from sources in terior to Neptune (Gomes 2003| 
|Morbidelli et al.||2003[ ). These models make the ad hoc 
assumption that hot and cold KBO populations have in- 
trinsically different colors (blue and red, respectively). 
This assumption is inconsistent with the data. The B — R 
distribution pays no regard to the reported hot/cold in- 
clination distribution. 

3.6. Double-checks 

Some of our objects possess large photometric error 
bars. Under the penalty of too low sampling we have 
also looked for the best fitting solution using only the 48 
CKBOs with colors having errors ^ 0.10. The best fitting 
solutions are found with equal probability for Eq.[7]when 
i k varies from 10.5° to 12.0° (with this sample we have 
no colors between these two inclination values). When 
comparing the mean colors of objects above and below 
i k with this reduced sample we also find equal incom- 
patibility levels for i c k — 10.5° up to 12.0°, as expected 
from the previous result. Also, when not discarding the 
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two apparent outliers mentioned in Section [2] the correla- 
tion values diminish slightly but all other results remain 
identical. 

Our sample of CKBOs was selected following a crite- 
rion based on the magnitude of the error bars instead 
of performing an average of all the published measure- 
ments for each object as used in the MBOSS sample (see 
Section [21 . To check the influence of this criterion we 
have double-checked our results using the CKBO colors 
from the MBOSS sample, since some color values were 
slightly different from those we have used. The outcome 
is identic al to that of Sec t ion 13.41 

Lastly, Gladman et al. ( 2008[ ) suggest an orbit classi- 
fica tion scheme slight!; 



different from the one used here 
The classifications of most 



(by iLykawka fc Mukai|2007 1 

KBOs remain unchanged between these two schemes 
and, not surprisingly, we find that our conclusions are 
statistically independent of the scheme employed. 

4. CONCLUSIONS 

The main goal of this work has been to investigate 
the color-inclination trend seen for Classical Kuiper Belt 
Objects (CKBOs) . We have analyzed a sample of B — R 
colors of 69 objects, excluding 2 apparent outliers as well 
as colors with error bars larger than 0.21. Obj ects were 
classified as CKB Os according to a definition by |Lykawka| 
& Mukai (20071. Orbital inclinations, denoted ik, were 
calculated relative to the Kuiper Belt Plane following 
Elliot et al. (2005). Our results may be summarized as: 



1. The linear B — R color-inclination correlation of 
CKBOs measured over the full range of inclinations 
from 0° to 34° is p = -0.70±g;{$, corresponding to 
a significance level larger than 8 a. This is a strong 



and highly significant correlation, consistent with 
previously published values. 

2. In contrast, the B — R colors of CKBOs with in- 



clinations ik ^ 12.0 



0+0.5 
-1.5 



are statistically uncorre- 
cted with inclination and are well described by 
B-R= 1.71 ±0.11. 

3. CKBOs with i k > 12.0 o t^ show a slight color 
vs. inclination dependence following (B — R) = 
—0.0159 ik + 1.703. The data are also formally con- 
sistent with a constant but bluer color, B — R = 
1.33 ±0.20, for ik ^ 12.5°, and a constant red color 
B-R = 1.70 ± 0.11 for i k < 12.5°. 

4. The data provide no evidence for a break or change 
in the B — R color distribution at the boundary 
between the dynamically hot and cold populations, 
purportedly near ik ~ 5°. In this sense, we find 
no observational support for the frequently-cited 
parity between red CKBOs and the dynamically 
cold population. The CKBOs are red up to ik ~ 
12° and, therefore equally red into the dynamically 
hot population. 
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Correlations and Wilcoxon Tests for Consecutive Inclination Cutoffs 
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a Orbital inclination cutoff Number of objects below (j) and above (h) the cutoff i£ c Mean B — R of objects 
below (i) and above (/,,) i£ B — R vs. if* correlation for objects below (;) and above (^) Significance level 
of the correlation^ Significance level of the Wilcoxon Test for color difference between objects below and above 
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TABLE 3 

Fit Results for Consecutive Inclination Cutoffs 



Two-constant stepwise fit b Constant-linear stepwise fit c 



il [°] a 


(B - R) ol 


(B - R) o2 


x 2 


(B - R) ol 


m 


(B - R) o2 


x 2 


3.0 


1.721 


1.528 


2.892 


1.721 


-0.0194 


1.801 


1.197 


3.5 


1.717 


1.508 


2.737 


1.717 


-0.0201 


1.817 


1.188 


4.0 


1.713 


1.489 


2.597 


1.713 


-0.0210 


1.837 


1.174 


4.5 


1.717 


1.473 


2.427 


1.717 


-0.0207 


1.830 


1.180 


5.0 


1.718 


1.465 


2.344 


1.718 


-0.0206 


1.827 


1.181 


5.5 


1.718 


1.465 


2.344 


1.718 


-0.0206 


1.827 


1.181 


6.0 


1.716 


1.459 


2.314 


1.716 


-0.0208 


1.834 


1.182 


6.5 


1.715 


1.444 


2.191 


1.715 


-0.0208 


1.834 


1.183 


7.0 


1.717 


1.434 


2.090 


1.717 


-0.0205 


1.825 


1.180 


7.5 


1.718 


1.423 


1.971 


1.718 


-0.0198 


1.809 


1.173 


8.0 


1.715 


1.417 


1.964 


1.715 


-0.0203 


1.821 


1.185 


8.5 


1.713 


1.410 


1.924 


1.713 


-0.0205 


1.825 


1.192 


9.0 


1.713 


1.410 


1.924 


1.713 


-0.0205 


1.825 


1.192 


9.5 


1.707 


1.395 


1.876 


1.707 


-0.0212 


1.844 


1.213 


10.0 


1.707 


1.395 


1.876 


1.707 


-0.0212 


1.844 


1.213 


10.5 


1.712 


1.374 


1.627 


1.712 


-0.0187 


1.777 


1.155 


11.0 


1.712 


1.374 


1.627 


1.712 


-0.0187 


1.777 


1.155 


11.5 


1.712 


1.374 


1.627 


1.712 


-0.0187 


1.777 


1.155 


12.0 


1.712 


1.345 


1.389 


1.712 


-0.0159 


1.703 


1.100 


12.5 


1.704 


1.328 


1.385 


1.704 


-0.0155 


1.691 


1.166 


13.0 


1.701 


1.316 


1.349 


1.701 


-0.0145 


1.665 


1.181 


13.5 


1.692 


1.319 


1.537 


1.692 


-0.0182 


1.765 


1.315 


14.0 


1.692 


1.319 


1.537 


1.692 


-0.0182 


1.765 


1.315 


14.5 


1.686 


1.296 


1.512 


1.686 


-0.0176 


1.748 


1.376 


15.0 


1.686 


1.296 


1.512 


1.686 


-0.0176 


1.748 


1.376 


15.5 


1.686 


1.296 


1.512 


1.686 


-0.0176 


1.748 


1.376 


16.0 


1.686 


1.296 


1.512 


1.686 


-0.0176 


1.748 


1.376 


16.5 


1.686 


1.296 


1.512 


1.686 


-0.0176 


1.748 


1.376 


17.0 


1.682 


1.282 


1.486 


1.682 


-0.0158 


1.696 


1.397 


17.5 


1.680 


1.265 


1.427 


1.680 


-0.0114 


1.570 


1.391 


18.0 


1.672 


1.265 


1.600 


1.672 


-0.0175 


1.746 


1.544 


18.5 


1.672 


1.265 


1.600 


1.672 


-0.0175 


1.746 


1.544 


19.0 


1.672 


1.265 


1.600 


1.672 


-0.0175 


1.746 


1.544 


19.5 


1.672 


1.265 


1.600 


1.672 


-0.0175 


1.746 


1.544 


20.0 


1.672 


1.265 


1.600 


1.672 


-0.0175 


1.746 


1.544 



a Critical orbital inclination value, i.e. location of stepwise behavior. Sec g3.4.2| and Eq.|4j c Sec 
|3.4.3| and Eq.[6] 



